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1 Introduction 



One of the simplest examples of a general problem of a quantum theory in a given external 
field is the study of quantum effects in the presence of mirror-like boundaries. A mirror 
can be considered as a (finite or infinite) potential barrier, the position of which is sharply 
localized in space. If mirror-like boundaries are static, the ground state of the system is 
still a vacuum, but the energy of this state is different from the energy of the vacuum in the 
absence of boundaries. This is the well-known Casimir effect A new phenomenon, 

the dynamical Casimir effect [0-0, occurs for moving boundaries. In the general case, 
a motion of mirror-like boundaries creates energy fluxes and particles from the vacuum. 
Non-uniformly moving mirrors radiate "photon" pairs. In many aspects, this effect is 
similar to the effect of parametric excitation of a quantum oscillator. 

Since the pioneering work of Moore |4]] who studied the creation of photons by moving 
mirrors for a 1-dimensional cavity, numerous publications investigating different aspects 
of the dynamical Casimir effect have appeared. Fulling and Davies found that the back 
reaction force generated by emission of photons in 2-dimensional spacetime (that is by a 
1-dimensional mirror) is proportional to the third derivative of the mirror displacement. 
Two-dimensional problems for conformal-invariant fields allow quite a complete study 
because of the presence of an infinite-dimensional group of conformal transformations. 
One can use these transformations to map the original problem with dynamical boundaries 
to a corresponding static problem which allows a solution. 

Absence of such a wide group of invariance makes the study of problems in higher 
dimensions much more complicated. Simple case when plane boundaries are moving with 
constant velocities was considered in || [|. Perturbation methods were used for study of 
physically interesting case of the emission of photons by a vibrating cavity in a realistic 
(3 + l)-dimensional spacetime JlO|-[0. This study demonstrates that the mirror-induced 
radiation, which is usually far too weak to be observed for speeds much below the speed 
of light, can be enhanced by orders of magnitudes by resonance effect due to the cavity 
finesse. A review of results for radiation from moving mirrors and bodies with finite 
refractive index can be found in [|18|]-[^(J . 



Not so much work has been done in the study of radiation created by a relativistically 
moving mirror in four (and higher) dimensional spacetimes. Candelas and Deutsch [^TJ 
considered radiation from a uniformly accelerated ideally-reflecting plane. The method 
of images was used in |E| and |23| to obtain the exact Green function and to calculate 
the stress-energy produced by a spherical ideal mirror, which is expanding in the four- 
dimensional spacetime with a constant acceleration for a scalar and electromagnetic fields, 
respectively. The latter calculations can be interesting in connection with the study of 
quantum effects produced by relativistic bubbles generated in the cosmological phase 
transitions (see e.g. - |25[|). 



In this paper, we continue the study of quantum effects produced by spherical mirrors 



2 



expanding with a uniform acceleration. More concretely, we study a scalar massless field 
propagating in a spacetime with partially transparent mirrors]]. The latter is described by 
5-like potentials. We consider two problems: Problem A for which the mirror potential 
is localized at the hyperboloid R 2 — T 2 = a~ 2 , where a is the value of the acceleration 
and R is the radius of the mirror at time T; and Problem B, where there exist two 
expanding^ concentric mirrors. This allows us to demonstrate how partial transparency 
of mirrors (which is a common feature of more realistic models) modifies the expression 
for the radiation emitted by expanding mirrors. After performing a Wick's rotation 
T — > iX , the problem is reduced to finding a potential generated by a point-like source 
in 4-dimensional space in the presence of one or two concentric 5-like spherical potentials. 
The latter problem can be solved by decomposing into spherical harmonics. Since the 
solution of the problem can be easily obtained in an arbitrary number of dimensions, we 
perform the calculations for the general case of Z)-dimensional spacetime. 

The paper is organized as follows. Section 2 discusses the formulation of the problem 
and considers two special models: Model A with one expanding partially transparent 
mirror, and Model B with two expanding concentric mirrors. Section 3 contains the 
calculation of the Euclidean Green functions for both problems. The expressions for (( 1 3 2 ) ren 
and (T^) ron for the Euclidean formulation are obtained in Sections 4 and 5, respectively. 
Radiation from a spherical mirror, which expands with a constant acceleration and a 
couple of expanding concentric mirrors, is calculated in Section 6. Section 7 contains a 
discussion of the obtained results. 

2 Formulation of the Problem 

2.1 General formulas 

Our aim is to study the influence of an accelerated motion of semi-transparent mirrors 
on quantum fields. In our consideration, the background geometry is flat. Nevertheless, 
it is convenient to first consider the general action for a scalar massless field in curved 
geometry 

S[ V ] = -\J dx D g 1 ' 2 (iTtpM + ZRtf?). (2.1) 

We keep the number of spacetime dimensions arbitrary and denote it by D. We also 
denote by d = D — 1 the number of spatial dimensions. The field equation is 

n<p-£R(p = 0, (2.2) 

1 For discussion of the Casimir effect for a system with partially transparent boundaries see, e.g. recent 
papers |^6|, and references therein. 

2 Strictly speaking for the motion with a constant acceleration, a spherical mirror at first contracts 
from an infinite size to a finite radius, where its velocity vanishes and only after this it starts its expansion, 
so that the velocity of this expansion reaches its asymptotic value c. 
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where □ = g l ^ 2 dAg 1 ^ 2 g^ u d u ). The stress-energy tensor, TL,, for the field has the form 
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T, v = (1 - 20 ^ ip, v + (2£ - - J ^ ^ ^ ^ 

- 2 ^ (</v„ - g^u <^;a v a ) + € (r^v - R) v 2 - (2.3) 

When the parameter £ of non-minimal coupling vanishes, reduces to the canonical 
stress-energy tensor. The case 

£ = -^ = ^ (24) 

* 4 D - 1 - 2d [ ' 

corresponds to the conformally invariant theory with = 0. 

Consider now the quantum field (p in a flat spacetime and denote by G (x,x') its 
Hadamard function 

G 1 (x, x') = (<p(x) <p(x') + <p(x') <p(x)) . (2.5) 
Then for a quantum average of stress-energy tensor we have 



(T^{x)) = \imD^G^{x,x'), (2.6) 

x — >x 

where 

1 /. 1\ ^ 



= - £J V M V,, + (£ - -J ^ <T V a 

- ^(V M V, + V^V,0- (2.7) 

Here V M and V M ' are gradient operators acting on the first and second arguments of 
G(x, x'), respectively. 

We write the metric of the flat spacetime metric in Cartesian coordinates as 

d 

ds 2 = -dT 2 + Y / (dX i ) 2 . (2.8) 

i=i 

We shall also use the Euclidean metric obtained from ( |2.8| ) by a Wick's rotation T — > iXq 

ds 2 E = j2(dX,) 2 . (2.9) 

The Green function G(x, x') of the quantum field ip in the flat spacetime is defined as 
the solution of the equation 

aG(x,x') = -6{x,x'), (2.10) 

obeying given boundary conditions. In the absence of boundaries, the Euclidean Green 
function Go is singled out by the condition that it vanishes at infinity and is written in 
the Cartesian coordinates (12. 91) as 



r f- 1 1 

G ^ *') = ^km 1 \x -xr- 2 ' (2 - U) 
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where 

\X-X'\* = J2(X,-Xlf. (2.12) 



This Green function is a solution of equation ( |2.10| ) with 



d Q2 

n * = E^- (2-13) 

The renormalized value of the stress-energy tensor is defined by ( |2.6| ), where instead 
of G^(x, x'), we must substitute 

G<M{x,x') = G^(x,x')-2G (x J x'). (2.14) 

Similarly, for (0 2 ) describing fluctuations of the field tp, we have 

(^(x)r = ilimG«(x ja 0. (2.15) 

Z 3T-+X 

2.2 Models 



To describe a semi-transparent mirror, we include in the equation ( J2.1U ) an additional 
term, — Vs 0, with a potential Vs of the following form: 

y s = y 5(S), (2.16) 

where E is a timelike d-dimensional hypersurface representing the motion of a (d — 1)- 
dimensional mirror surface. For finite values of Vq, a mirror can be partially penetrated 
by the field ip. In the limit Vq — > oo the penetration coefficient vanishes, and one has an 
ideal mirror which reflects all the frequencies. 

We consider two models of moving semi-transparent mirrors. 

2.2.1 Model A 

As a first example, we consider the case when there exists a single mirror and its surface 
worldsheet E is described by the equation 

E : R 2 -T 2 = b 2 , (2.17) 

where 

R 2 = J2xf. (2.18) 

i=l 

Such a mirror is spherical. Its radius R changes with time T in such a way that points of 
the mirror are moving with a constant acceleration a = b^ 1 orthogonal to the surface of 
the mirror. In (T, i?)-coordinates the motion is represented by a hyperboloid (see Fig.l). 
The Wick's rotation T — > iXq transforms E into E# 

E E : R 2 + X 2 = b 2 , (2.19) 

which is a <i-dimensional sphere S d of radius b embedded into D-dimensional Euclidean 
space. 
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Figure 1: The geometry of problem A. Solid lines in the left figure represent motion of spherical 
semi-transparent mirror £ in the spacetime. A dash-dotted line in the right figure corresponds 
to the surface Ee in the Euclidean space obtained by the Wick's rotation T — ► iXq. 




Figure 2: The geometry of problem B. Solid lines in the left figure represent motion of spherical 
semi-transparent mirrors Si and £2 m the spacetime. Dash-dotted lines in the right figure 
correspond to the surfaces Si._e and ^2,e in the Euclidean space obtained by the Wick's rotation 
T - iX . 

2.2.2 Model B 

Another example is a case when there exist a two concentric spherical mirrors. In this 
model, the surface S consists of two concentric mirror surfaces, Si and E 2 , given by (see 
Fig.2) 

51 : R 2 -T 2 = bl, (2.20) 

5 2 : R 2 -T 2 = b\. (2.21) 

We assume that 62 > b\ and call S2 and Si an external and an internal mirror, respectively. 
After the Wick's rotation T — > iX , these surfaces become (see Fig.2) 

51, £ : R 2 + X 2 = b\ , (2.22) 

5 2 , £ : R 2 + X 2 = b 2 2 . (2.23) 
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3 Calculation of Green Functions 



3.1 Green Functions in Spherical Coordinates 

We start with the Euclidean formulation. It is convenient to introduce spherical coordi- 
nates in the ^-dimensional^ Euclidean space (|2.9|) . For this purpose we first parameterize 
the line element dQ d on a unit d-dimensional sphere as follows: 

dnl = d$l, (3.1) 

dtt 2 n = d&l + sin 2 d n dtt 2 ^ , (n = 2, . . . , d) . (3.2) 

In these co-ordinates 

d 

dtt 2 = Y, n f d&i dtfj , (3.3) 

n[f =diag(l,sin 2 i? d (l,sin 2 i? d _ 1 (... .))) , (3.4) 



/n d = ^det [plf) = sin^" 1 $ d sin^ 2 . . . sin tf 2 . (3.5) 
We can write the line element (12.91) as 



ds\ = dp 2 + p 2 dQ 2 d , (3.6) 

where 

R = psm-dd, Xq = p cos $d . (3.7) 
In what follows, we shall also use the notation 

V = #d. (3-8) 
In spherical coordinates, the Euclidean "Box" -operator D E takes the form 



d ^|(4) + >' (3 - 9) 



where A d is the Laplace-Beltrami operator on S d . 

Scalar spherical harmonics on S d , which are solutions of the eigenvalue problem 

A d Y = XY (3.10) 



on S , have been studied by many authors (see, e.g. |28|, |29|, |30|). Let li,l%, • • • , Id be 
integers that satisfy 

Id > k-i > ••• > k > \h\, (3.11) 



3 We assume that D > 2. A special case of a two-dimensional spacetime which requires slight 
modifications will be discussed in Section 5.3.3. 
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and denote 



2\ + n-l (X + l + n-2)\ 



1/2 



(sin??) 



-(n-2)/2 



(A-/)! 

A+(n-2)/2 V-OS^), (3.12) 

where P£(x) is the associated Legendre function of the first kind. Then the spherical 
harmonics on S d are given by 



Ln=2 



'27T 



(3.13) 



Scalar spherical harmonics are normalized as 

d#! . . . d# d V^S V- ii ' %- *i = ^ i' d • • • *Wi > ( 3 - 14 ) 

and form a complete set of functions on 5* d . Scalar spherical harmonics obey the relation 

& d Yi d ... h = -l d {l d + d-l)Y ld ... h . (3.15) 
To make the notations brief, we denote 

L = l d , W = {l d - 1 ,...,l 1 }> (3.16) 



where values of l n , which enter the collective index W, are restricted by ( |3 . 1 1 ) . In these 
notations equations ( BTT5 ) and ( |3.14| ) take the form 



Here 



and 



A d Y LW = —L(L + d- 1)Y LW , (L > 0), 
J dVL d Y LW {yt) Y L > W i{p!) = 5 L l' 5ww> ■ 

dVt d = d& 1 ... d'd d \j£l d , 
$ww = K-ii' d ^ ■ ■ ■ Ki[ ■ 



We also have the relation 



Y,Y, Y Lw(n)Y LW (n') = 5(n,n>) 

L=0 W 



(3.17) 
(3.18) 

(3.19) 

(3.20) 

(3.21) 



which is a consequence of the completeness of scalar spherical harmonics, and where 
5(p., Q f ) is an invariant 5-function on unit sphere S . 

The external potential Vs B for both of the problems, A and B, is evidently invariant 
under rotations in the Euclidean space. That is why both problems allow for separation 
of variables. Substituting (|3.21|) into the equation 

(n E - Vx E )G E (x, x') = -5{x, x') , (3.22) 
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we can obtain the following representation for the Euclidean Green function: 



g e (x,x') = ^T,^L(p,p')Y LW (n)Y LW (n'). 

L=0 W 



(3.23) 



The radial Green function Gl(p,p') obeys the equation 



d ( d d 



p^L(L + d-l)-p d Vx(p) 



g L (p,fO = -8(p-ft)- (3-24) 



dp \ dp j 

Equation (|3.24j) is self-adjoint, and hence Gl is a symmetric function of its arguments 

Gl(p,p') = Gl(p',p). (3.25) 

For problem A 



while for problem B 



Vx(p) = V 6(p-b) 



V s {p) = V l 8{p-b 1 ) + V 2 5{p-b 2 ). 



(3.26) 



(3.27) 



The radial Green function Gl does not depend on the collective index W. One can use 
the following relation to make the summation over W in (|3.23|) ( ||29|| , equation (B.12)): 



2 Y LW {Q) Y LW (Q') = r]d (2L + d- 1) ci^ 1)/2 (cos 7 ) 



(3.28) 



where 7 is the angle between fl and fl' on the unit sphere S , and 



Vd 



47r (d+l)/2 • 

Note that the volume of the d-dimensional unit sphere is 



(3.29) 



and hence 



The function cf 1 I 
function F(a, b; c; z) as follows 



2tt 2 



r(¥) 



(d-l)V d 



(3.30) 



(3.31; 



x) is the Gegenbauer polynomial related to the hypergeometric 



(L + d-2)\ 



d 1 — x s 



C ( L d - 1)/2 (x) = ^—f — =f F(-L, L + d- 1;-; 
L v ; L!(d-2)! v ' 2' 2 



(3.32) 



Since the hypergeometric function is normalized so that F(a, b; c; z = 0) = 1, we get 



C 



(d-l)/2. 



(L + d- 2)! 
L\ {d-2)\ 



^3.33) 
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Thus, we have as the following representation for Ge- 



G E {x, x') = Vd J2(2L + d- l)g L (p, p') Cr i)/2 (cos 7 ) • 



(3.34) 



L=0 



In this expression, 7 is a geodesic distance between two points . . . d d ) and (i?-^ . . . i? d ) 
on the unit sphere. It is defined as 7 = 7^ by the following relations 

71 = A - #i , (3-35) 

cos 7„ = cos $ n cos d n + sin "d n sin $ n cos 7 n _i , n = 2, . . . , d . (3.36) 

An ambiguity in the choice of the spherical coordinates connected with rigid rotations of 
space can be used to put $1 = . . . = "Qd-i = and ■&[ = ... = $ d _i = 0) for any chosen 
pair of points on the unit sphere. 

3.2 Problem A 



In order to construct Ql{p,p') for problem A , we shall use solutions of the homogeneous 
equation 

\ J ( . J \ 

R L = 0. (3.37) 



U" d h)-^ L{L+d - i: 



It is easy to verify that this equation has the two linearly independent solutions, p L+ and 
p L - , where 

L + = L, L_ = -(L + d-1). (3.38) 

The solution p L+ is regular at the origin p = 0, while the solution p L ~ decreases to zero 
at infinity. The Green function Gl(p,p') which is regular at p = and decreases to zero 
at infinity is of the form 



' (A- p L - + A + p L+ )(p') L+ , b> p> p' >0, 

Bp L - (p') L + , p > b > p' > , 

4 P L - (C_ ( P ') L - + C + (p') L+ ) , P > P' > b > . 



(3.39) 



The symmetry condition ( 3.25|) determines Ql{p,p') for other possible relative positions 
of p and p' . 

The radial Green function Gl{p, p') is continuous, while its first derivatives have jumps. 
The form and location of the jumps directly follow from equation (|3.24 ) with potential 
( 13.261 ). The jump at p = p' is of the form 



dg L {p,p') 




dg L (p,p') 


dp 


p=p'+0 


dp 



p=p'-0 



(p 



(3.40) 
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For p' ^ 6, the jump at p = b is 



dGL(p,p') 




dGL(p,p') 




dp 


p=b+0 


dp 


p=b~0 



v g L (b,f/). 



(3.41] 



The jump at p' = b and p ^ b has a similar form. Relation (|3.40|) allows us to find 
coefficients A_ and C + in ( |3.39| ). Continuity of Gl(p, p 1 ) aX p = b and p' = b, together with 
relation (|3.41 ), fix the other coefficients. By solving the corresponding linear equations, 
we get 

A_ = C + - ., .. " . . D 



1 



CL 



2(L + (3 d ) 



2(L + f3 d + U ) 
1 



2(L + /3 d )(L + (3 d + U )b^+d-i 
U 



(3.42) 



2(L + f3 d )(L + (3 d + U ) 



h 2L+ 



d-l 



Here we introduced notations 



fa = d -^- , U = l -bv Q 



(3.43) 



In the absence of a mirror, when Vq = 0, relations (|3.39|) are simplified to become 

1 



A_ = C , = B 



2L + d-l 



, A + = C_ 



0. 



(3.44) 



In this limit, the Green function Ge(x } x') coincides with the free Green function (|2.11 
This can be easily verified by using the relation (see e.g. p9| , eq. (B.13)) 

1 



\X -X 



nd - 1 Vd^ ^2L + d-l 



(3.45) 



Here p > = max (p, p') and p< = min (p, p'). 

Subtracting the free Green function Gq(x, x') from Ge(x,x') gives the renormalized 
Green function 



GT(x, x') = y, E sr (p, P') YlwW Y lw (SY) , 
l=o w 

= % £(2L + d-l) Grip, P') ^r i)/2 (cos( 7 )) . 



(3.46) 



Here 



Gr(pip') = - 



Ua 



2(L + (3 d )(L + (3 b + U ) ¥-i 



EL 

V 
6 2 



, b > p, p> > 



P, P'>b- 



(3.47) 
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3.3 Problem B 



The calculation of the Green function for the problem B is similar to the calculations 
of the previous subsection. We again use the representation ( |3.23| ) and write the radial 
Green function Gl(p,p') m the form[] 



' (a_ p l - + A + /+) (p'Y 



B_ p L - +B +P L +)(p') L + 



Cp L -(p') L + 



Gl(p,p') 



h > p > p' > , 

b 2 > p > h > p' > 0, 
p > b 2 > bi > p' > 



C- 1 (p_ p L - + D + p L +) (E4p') L - + E+(p') L +), b 2 > p > pf > h > 



P- (F4p') L - + F + (p') L+ ) 
( p l - {G-(p') L -+G + (pr 



p > b 2 > p' > bi > 



p > p' > b 2 > h > . 

(3.48) 

The expression for Ql(p,p' for other possible positions of p and p' follows from the sym- 
metry of this function. As before, jump conditions ( |3.4ip and ( |3.42| ) allow us to obtain 
the unknown coefficients which enter relation (|3.48|) . After straightforward but long cal- 
culations we find that 



1 



2Q L bj L+d ~ 1 
A_ = G+- 



1 1 --^-(1-7/J 



1 



B+ 


= D+ 




= F+ 



2(L + &)' 
U 2 1 



2n L bf^- 1 ' 



G- 



Here, we used notations 



2fl L 

h 2L+d-l 
^2 

2ttr, 



b 2 V 2 



C = 
B_ = D_ 



L + (3 d 



L + (3 d + U 2 



2Q 



2n L 
-K 1 



Ul i2L+d-l 



U llL + U 2 +^-(l- lL ) 



1L 



2L+d-l 



Pd 



d- 1 



(3.49) 

(3.50) 
(3.51) 
(3.52) 
(3.53) 

(3.54) 
(3.55) 



n L = (L + [3d + U x ) (L + (3 d + U 2 ) - U x U 2lL . 

For calculation of (0 2 ) Icn and (T^) ron we need £7™, which can be obtained from (|3.48| ) 
by subtracting the vacuum radial Green function Q\. The latter can be obtained by 

4 Coefficients A±, B±, . . . which enter relations ( 3.48 ) - ( 3.53 ), ( |3.58| ) and ( 3.59 ), depend on the 
multipole moment L. In order to simplify formulas we do not indicate this dependence explicitly. 
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setting U\ — U 2 = in ( |3.48| ). In the regions of interest (p' < p < b\ (inner), p > p' > b 2 
(outer), and b 2 > p > p' > b x (intermediate)) the vacuum radial Green function has the 
same form 

i 

(3.56) 



1 



2(L + /3 d) 

As a result, we get the following expressions for Q T £ n : 



L_ iLj, 
p p 



Inner region p,p' < b\\ 
Outer region p, p' > b 2 . 



l\L 



ST(p,p' 



G. 



(pp 



/\L+d-l 



(3.57) 



(3.58) 



Intermediate region b 2 > p, p' > b\\ 

1 



2il L (L + f3 d 



2L+d-l 



(pp 



+U 2 (L + (3 d + U x ) -^zt " U! U 2 lL 



(P 



^pi^L+d-l pL+d-1 



(3.59) 



As expected, {?£ cn (p, p') is a symmetric function of its arguments. It is easy to verify 
that when one of the potentials, U\ or U 2 , vanishes, expressions ( |3.571 )- (|3.59| ) for Q r £ n 
reduce to (|3.47|) . In the other limit of perfectly reflecting mirrors, when U\ = U 2 = oo, 
relations ( |3.57| )-( p.59p correctly reproduce the result of . 



4 Calculation of ((p 2 ) ren 

4.1 Problem A 

4.1.1 Expression for ((p 2 ) ren 

In order to calculate (0 2 ) Icn , we need to obtain G^ n (x, x') in the coincidence limit x' — > x 

m- = GT(x,x). (4.1) 



Let us discuss problem A at first. We shall use representation (|3.34|) for the Green 
function, where, instead of Ql, we substitute its renormalized version Q T £ n given by (|3.47|) . 
In order to find a coincidence limit, we must put 7 = and p = p'. Using relations (|3.33| ), 
( PD , and we get 



mp)r 



UoVd 
b d-i 



FW((p/b) 2 ,U + (d- l)/2), p<b. 
(ff {d ' 1] F^((b/ P r,U + (d-l)/2). b>p; 



(4.2) 
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where the function is 



F«(z,/3) = X; 



(L + d-2)! 



(4.3) 



^ L! (d-2)!(L + 0) 

Some properties of this function are discussed in Appendix A. In particular, relation (|A.3| ) 
allows one to show that at the center p = 

U 



(<p'(p = o)p = - 



Vd 



(4.4) 



b^ 1 (d - l)/2 + U ' 

Quantity (<^ 2 (p)) rcn is divergent at p = b. Using ( |A.5| ) one gets the following expression 
for the leading divergent part of (0 2 {p)) Icn at p = b 

U Tj d 1 



mp)Y 



where 



W- 1 (d - 2) (1 - x 2 ) d - 2 
p/b, p<b, 



x 



(4.5) 



(4.6) 



b/p, b> p. 

In the limiting case Uq = oo, (ip 2 (x)) Icn can be expressed in terms of elementary 
functions. Using ( |A.6[ ) we can rewrite ( [4.2| ) for Uq = oo in the form 

Vd 1 



b d - x |1 - ( p /6) 



2 Id— 1 



(4.7) 



For d = 3, this result directly follows from expression (|3.5| ) of Ref. |22| for the Green 
function of a scalar massless field in the presence of an ideal mirror with the same choice 
of the surface £ as for our problem A. For d ^ 3 the result (|4.7|) can be easily verified by 
the method of images which was used in Ref. ||22|| . 

4.1.2 4-dimensional case 

Now we discuss some properties of (0 2 ) rcn in the case of the 4-dimensional spacetime, 
which has the most physical interest. Figure 3 shows plots of (0 2 ) Icn (for d = 3) as a 
function of p/b for different values of the height of the potential Uq. For comparison, it 
also contains (<p 2 ) rcn for an ideally reflecting mirror, (|4.7Q . 

It is evident that, near the mirror surface, the dominant divergent part is directly 
related to the singularity of the potential. In order to extract the singular part of (<p 2 ) rcn , 
we represent in the form 



F^(z,f3) = f(z,(3)+g(z,(3) 



(4.8) 



where 



h(z,/3) 



1 - z 



;i-/3)ln(l-z) + L^ + 



1+/3 



"(2-/?) 



z + 



2 + p 



+ 



(3-3 



z . 
(4.9) 
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0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 



Figure 3: — 4ir 2 b 2 ((p 2 ) lan in 4-dimcnsional spacetime as the function of p/b for different values 
of the parameter Uq: solid line - Uq = 0.2, dotted line - Uq — 1.0, and dashed line - Uq = 5.0. 
The dashed and dotted line corresponds to an ideally reflecting mirror (Uq = oo). 




Figure 4: Function H(p,Uo) for different values of the parameter Uq: solid line - Uq = 0.2, 
dotted line - Uq = 1.0, and dashed line - Uq — 5.0. 



1.4- 




Figure 5: Function G(p,Uo) for different values of the parameter Uq: solid line - Uq — 0.2, 
dotted line - Uq = 1.0, and dashed line - Uq = 5.0. 
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g(z,f3) = -(l-(3)~g(z,f3) 



{L(L + (3) l-(3 



l + (3 



(2-/3) 



z + 



3 (3-3 
+ 



2 + (3 



z 2 . (4.10) 



The series in the right-hand side of ( |4.10| ) is convergent at z = 1. Hence, the function 
f(z, (3) contains all the divergences of at z — 1. The last terms in the right-hand side 
of ( |4.9| ) which are proportional to z and z 2 are finite at z — 1. They are added in order 
for /i(z, /3) to have the same first 3 terms in the Taylor series expansion at z = as the 
original function F^ 3 \ Hence, we can write (<£> 2 ) rcn in the form 



2U2 



where 



4n 2 b 



H(p,U ) 



G(p,U ) 



■&(p) 2 r n = H(p,U ) + U G(p,U 



h((p/b) 2 A + U ), 

(b/ P n((b/p) 2 , 1 + U ) 
~g((p/b) 2 ,l + U ), 



P<b, 

b> p; 
P<b, 



(4.11) 



(4.12) 



(4.13) 



{b/pf~g{{b/pf,l + U ), b>p. 
The plots of functions H and G for different values of U are shown in Figures f| and ^[ 

4.2 Problem B 



Using relations (|4.1| ), ( |3.33| ), and (|3.57|) - (|3.59|) , we can write the following representation 
for (ip 2 ) rcn for the two-mirror case (problem B): 



Vd g (L + (3 d )(L + d-2)\ 



(d-2 



L=0 



(4.14) 



where 



K L (p) 



1 f P_\ 2L 
^Lbf- 1 \bj 



n L (L + f3 d ) P d ~ l 



n -l r -i-L (1"7l) 



P < b 



i ■ 



V 



2L+d-l 



(4.15) 



+U 2 {L + p d + U 1 ) 



\ 2L+d-l 



2L+d-l r 



and 7L = {b 1 /b 2 ) 2L+d - 1 . 



U 2 + U llL + 



2U x U 2lL 

U1U2 ( 

L + f3 d 1 



1-7l) 



b\< p <b 2 . 
, P > b 2 . 
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bl=1.0, b2=3.0, Ul=100, U2=100 




Figure 6: — Air 2 b 2 {(p 2 ) TCU in 4-dimensional spacetime as the function of p/b for different values 
of the parameter U\ and \Ji. 

Notice that only the L = modes contribute to (<^ 2 ) ron at the origin. Thus, one has 

where 70 = (bi/b 2 ) d ~ 1 ■ 

To illustrate a typical behavior of (y3 2 ) rcn for the two-mirror problem, we plot in Fig- 
ure 6 the value of — 4ir 2 b 2 {(p 2 ) icn in 4-dimensional spacetime as a function of p for different 
values of the parameter U\ and U2, and for b\ = 1 and 62 = 3. 



5 Stress-Energy Tensor 
5.1 General properties 

For both problems A and B, the boundary conditions are invariant under the group 
of rotations 0(d + 1) which is the space symmetry. As a result of this symmetry, the 
renormalized stress-energy tensor has the structure 

(f;) rcn = diag(e,p,...p), (5.1) 

where e and p are functions of p. The conservation law 

{f;T% = (5.2) 
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implies 



de d 

— = -- (e-p . 

dp p 



(5.3) 



For the conformal invariant theory (|2.3|) — Q2.4|) we have (X^) rcn = e + d ■ p = 0, while 
the conservation law (|5.3j ) gives 



_eo_ 



p = - . (5.4) 

In other words, the renormalized stress-energy tensor for the conformal invariant theory 
is uniquely determined by the symmetry and the conservation law up to one constant. In 
a more general case, it is sufficient to determine only one function of one variable p, say 
e — p. 

Using equations ( |2.6| ) and (|2.7|) , we get 



e — p = lim 



p 2 



G rcn {x,x'). 



Here 



K = (1 - 2£)d p d p , - £ (d p 2 + d 2 p ) +- p (d P + d pl ) 



(5.5) 
(5.6) 

M = - (1 - 20 <9„ + £ [d 2 v + cfy] . (5.7) 

(The last term in the right-hand-side of (|5.6|) arises because in the spherical coordinates 
I^L = — p 7^ 0.) We use now this relation to obtain the renormalized stress-energy tensors 
for our problems. 



5.2 Problem A 

Using expressions ( |3.46| ), ( |3.47[ ), and (p . 5|) , we can write 

U T]d 



e — p 



7 = 



where 



Q±(ft "' 7)= £(rw^) ^ 



COS 7) . 



It is easy to verify that 

\n(p P ') L± 



p=p' 



;i-40 L 2 ± + 4£L ± 



P 



2(L±-1) 



A/-cl d - 1)/2 , 



cos 7) 



C 



(d-l)/2. 
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(5.8) 



(5.9) 



(5.10) 



(5.1i; 



7=0 d 

In order to obtain ( |5.11| ), we used ( |3.32| ) and the following property of hypergeometric 
functions 



d 

— F(a,b; c; z) 



2=0 



ab 
c 



(5.12) 
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Using relations ( |5.10D , ( |5.11| ), and ( |3.33j ), we get 



UoVd ^ 



I (L + d-2)\ (p 2 \ 



2 \ L ± 



W-ip* ^ (L + (3 d + U Q ) L\{d-2)\ \V) 



(5.13) 



where 



(5.14) 



The signs + and — correspond to the internal (p < b) and external (p > b) problems, 
respectively. Using relation ( |5.3| ), we obtain the following expression for e 



UoVdd ~ Aj (L + d-2)\ (p 2 \ L± 

2¥- 1 p^ (L ± -l)(L + (3 d +U ) L\(d-2)\ \V ) 



(5.15) 



It is easy to see that first two terms (with L = and L — 1) ( 5.13j ) for the inner problem 
vanish. Hence, in the general case, e — p ~ p 2 for small p, while e remains finite at p = 0. 



Equation (|5.14j) shows that, for a special case of conformal theory, when £ = £d 



d-l 
Ad ' 



the coefficients A L ^ d vanish. For this reason, e = p, and hence e = p = 0. This result is 



similar to what happens in the case of an ideally reflecting mirror considered in p2 
When £ = 0, we have 



A^ =0 = ^ L ± (L ± - 1) . 
Let us consider the internal (p < b) problem first. In this case, 

(d-l)L(L-l) 

Substituting A\^ =0 m to (|5.13| ), we find for p < b 

e _ p =_ BUo jL F (^2) ( (p/6)2? ^ + {d + 3)/2 ) 

where is the function determined by relation (|4.3| ), and 

B =^(d-l) 2 . 

Similarly, for the external (p > b) problem, 

_ (d- 1)(L + d)(L + d - 1) 
" d ' 

Substituting this expression into fl5.13| ), we get 

b d-i 



(5.16) 



(5.17) 



(5.18) 



(5.19) 



(5.20) 



e — p 



BU 



id+2 



P 



2d 



(b/ P y, u Q + (d-i)/2 



(5.21) 



19 



Equation ( |5.3| ) allows us to obtain e (and hence p). Using relations ( |A.8| ) and ( |A.9| ) 

we get 

f F^ l \Up/b)\ U + (d + l)/2)), p<b, 



BUp 
2W+ 1 



J', 



Id 

| F^(((b/ P f,U Q +(d-l)/2)), p>b. 



(5.22) 



For an ideal reflecting mirror, Uq = oo, the expressions for e — p and e are greatly 
simplified. Using relation ( |A.7| ) we get 



kd-l 



e = ± 



B_ 

2 \b 2 -p 2 \ d1 

2 



\b 2 - p 2 



(5.23) 



(5.24) 



The signs ± in ( [5.23D correspond to the inner and outer problem, respectively. For d = 3, 
expressions Q5.23D - ( |5.24| ) reproduce the results obtained in |22| . 



5.3 Problem B 

5.3.1 General formulas 

Calculations of e— p and e for problem B are similar to the calculations for the problem A, 
but are much more involved. Besides relations ( |5.10| ) and ( |5.11| ), we use also the relations 



K (p L ± p' 1 



p=p p' 



4+1 



[L + L_-£(d 2 -l)]. 



For this section, we only list the final results. 
Inner region p < b±: 



Vd ~ (L + p d )(L + d-2)\ 
P ~ p> to L\{d-2)\ A ^'^P^ 



Vd d ~ (L + (3 d ) (L + d-2)\ 
' ^ to L\{d-2)\{L-1) A ^ M 



2p 2 



nip) = 

Outer region p > b 2 : 

e — p 



P_ 

bt l Ui 



2L 



U 1 + U 2lL + ^-{l- lL ) 
(L + p d + U X ){L + p d + U 2 ) - U,U 2lL 



rjabt 1 j^(_L + f3 d )(L + d-2)\ 



Q 2d 
r L=0 



L\ (d-2) 



Vdbt'd ^ (L + fo) (L + d-2)\ 



2p 2 * f^o L\(d-2)\(L + d) 



(5.25) 



(5.26) 
(5.27) 

(5.28) 



(5.29) 
(5.30) 
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2L 



Intermediate region b 2 > p > b\ 



f/i7L + f/ 2 + ^^-(l-7L) 
^ + Pd 

(L + f3 d + U X ){L + (3 d + U 2 ) - U x U 2lL 



(5.31) 



' P ~ P d+1 h L\(d-2)\ 



(5.32) 



_ Vd d ~ (L + d-2)! 
' Z= 2/+i ^ L!(d-2)! 



(L_ - i; 



.4 



f 40 7?o 



(5.33) 



The coefficients Af, in relations (|5T26| ) , ( gjg) and ( fxgp are given by equation (|5.14j) 



and 



^ = ^L + L_-^ 2 -l) = - 



d 



d+ 1 
d 



L(L + d- 1) + £(d 2 - 1) 



We also use the notations 

6z(p) 



U x {L + p d + U 2 



(L + ft + C/i)(L + ft + C/ 2 ) - U x U 2lL \ p) 



2L+d-l 



U 2 {L + (3 d + U 1 ) 



K 



P_ 

(L + ft + Ui)(L + ft + U 2 ) - U x U 2lL \h 2 



2L+d-l 



(5.34) 

(5.35) 

(5.36) 
(5.37) 



" L (L + ft + U 1 )(L + ft + U 2 ) - U x U 2lL ' 
Expressions for e in each of the three regions can be easily obtained by integrating relation 



5.3.2 Conformally invariant theory 



Equation ( |5. 14 ) implies that A^^ = for the conformal invariant theory, when £ = £ d 



(d - I) /Ad. Equations ( p37D , ggTgg ), ( gg^ ), and flOO]) show that e = p = both in the 
inner and outer regions. In the intermediate region, we have 

C 



e — p 



where 



d+1 « (L + d-2)! 
C =2UiU 2 r] d — 3— 2^ 



P' 



1L 



,d+l 



(5.38) 



L(L + d- 1) + 



(d-1) 



d ^ L!(d-2)! (L + ft + £/i)(L + ft + £/ 2 )- 
Using relation e = —dp we get 

e 



d C 



(5.39) 



(5.40) 



d + 1 p d+1 ' 

It is easy to verify that in the limit of ideally reflecting mirrors (U\,U 2 — > 00) in 4- 
dimensional space (d = 3), relations (|5.38|) - (|5.40| ) correctly reproduce the results of 



221. 
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5.3.3 Two-dimensional case 

We now make several remarks concerning (T^) ron calculations in a 2- dimensional space- 
timeQ. Since £d=% = 0, conformal and canonical stress-energy tensors are identical in 
2-dimensional spacetime. We restrict ourselves by discussing only this special case and 
put £ = 0. In the 2-dimensional case, the spherical harmonics (which enter for example 
into the expression ( |3.23| ) for the Green function) are simply functions exp(imrj)/ V2n, 



and (L, W) = m = 0, ±1, ±2, .... For all modes with m^O, the general relation ( |3.59| ) 
can be easily verified for the 2-dimensional case. (One need only to put (3d = there.) 
The only problem is due to the mode m = 0. A general solution of the radial equation 
( p. 24 ) for m = is logarithmically divergent either at p = or p = oo. For this reason, 



in the general case there does not exist a radial Green function Q m=0 (p, p') which remains 
finite at both boundaries. Fortunately, these zero modes does not contribute to (T^) rcn 
and the corresponding ambiguity is not important for our problem. The calculations for 
the problem B give e = p = inside the inner mirror and outside the outer mirror, while 
between the mirrors we have 



2U 1 U 2 ~ m 7m 

11 F m=l i6m 

C = (5 - 41) 

Here, 7 m = (&i/&2) 2m and Q m = (in + U\){m + f/ 2 ) — U\ Ui^i m ■ Both quantities e and 
p evidently vanish if one of the potentials Ui vanishes. This is exactly the result which 
must be expected for Problem A. 

6 Radiation From Mirrors 

Let us now discuss the properties of quantum radiation created by spherical mirrors which 
expand with a constant acceleration (see footnote 1). Our primary aim is to derive an 
expression for radiation at infinity produced by such a mirror. This problem can be solved 
by making an analytical continuation of the results of calculations for (0 2 ) Tcn and (T^) ren 
in the Euclidean space. 

6.1 (<£ 2 ) ren at J + 

First, let us study (0 2 ) rcn at J + . We consider only Problem B, since the result for problem 
A can be easily obtained by taking the limit U\ — ► 0. For Problem B, in the exterior 
region we have (cf. ( f4.14|) and (|5.31| )) 



5 We shall not consider (<^ 2 ) ren because of the infrared problem it is not a well-behaved quantity. In 
particular, (ip 2 ) rcn is logarithmically divergent at infinity for any value of the mirror potential. 
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Using (p.7|), we can write 

p 2 = R 2 + X 2 . (6.2) 
After performing a Wick's rotation X = iT, this quantity becomes 

p 2 = R 2 - T 2 . (6.3) 

By letting w = T — i? be a retarded time, we have at large R 

p 2 ^-2uR. (6.4) 

For this definition of u, the value u = corresponds to the moment of the retarded time 
when both expanding mirrors reach J + . For a part of J~ + lying to the past of this moment 
(that is outside the mirror) we have u < 0. The leading contribution at J + is given by 
the term L = in Q6.ll) . Hence we have 

(^)-(u) « - ^ r— 7^ , (6.5) 

\r / \ J R d-i (_2 M )«i-i 2 ' v ; 

0i7o + C^2 + ^ (1 - 7o) 

7?- — £s C« fil 

{Pd + U x ){p d + U 2 )-U x Uxy Q ' K ' } 



where 70 = (61/62) 



d-1 



6.2 (T^) ren at J+ 

Our starting point is the general expression for (T^) ron , which in the Euclidean space 
under consideration takes the form 

T> = <tT n = (e~p) 4 P) % + P% ■ ( 6 - 7 ) 

Here is the Kronecker symbol which is non-vanishing only when v is the index corre- 
sponding to the co-ordinate p. (No summation over p!) Using expression (|6.7D, we get in 
(X , R) coordinates 

Tx Xo = -5 (Xq6 + R 2 p s ) , 

Tx r = \x R{t-p), (6.8) 
P 



Trr = \(R 2 e + X 2 p) . 



P 

By making a Wick's rotation X = iT, we get 



T TT = ±( T 2 e-R 2 p), 



p« 

T TR = ~TR(e-p), (6.9) 
P 2 



Trr = - 2 (R 2 e-T 2 p 
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For fixed u = T — R and R — > oo we find 



R 2 

T TT « —Ttr ^T RR ^—(e-p). (6.10) 

P 



For large i?, the leading contribution to e — p is given by the L = term in (|5.29 ). Thus, 
we have 

p 2d 2 

where is given by (|6.6|) and 



t-P*-Hi ^o^o, (6-11) 



A^ = {d-1) (d-l-4fd) . (6.12) 

Combining these results, we get for T^ v in the asymptotic region the following expres- 
sion 

E 

T^ u = — ' ( 6 - 13 ) 
where = w iAt , u = T — R is retarded time, 

* S f - -I B (^ 5) (=S« ^ — l) 2 (d — 1 — 4 5 ,) n, , (6,4) 

and 

^d-^^Vd-i (6.15) 

is the surface area of a (d — 1) -dimensional sphere of radius R. The function B(z,w) in 
relation fl6.14|) is the Beta function B(z,w) = T(z)T(w)/T(z + w). 

For an ideally reflecting external mirror, U2 = 00, TZq = 2/(d — 1), so that for the 
radiation of such a mirror in 4-dimensional spacetime (d = 3) we have 

dE hi , , 

*r = -4^< 1 - 6 «' <" 6) 

For £ = 0, this result reproduces the result obtained in . 

By making the substitution u^v = T + Rm relations ( |6.5|) and (|6.13| )- ( |6.14|) we 
easily obtain expressions for (y3 2 ) ron and T^ v at J~ . 

In a similar way it is easy to obtain expressions for (0 2 ) Icn and T^ v in the region located 
inside the inner mirror or in the region between the mirrors (for problem B). We do not 
reproduce here the corresponding results. We only mention that both (0 2 ) rcn and T^ v 
remain finite on the null cone N : R 2 = T 2 for any finite value of time T. This property 
directly follows from the finiteness of (ip 2 ) ren and at the origin p = in the Euclidean 
formulation. The null cone iV plays the role of the event horizon for an observer moving 
on the mirror surface. Regularity of physical observables on the horizon differs the case 



of moving mirrors from the case of the black hole (see e.g. discussion in [[33 
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7 Conclusion 



In this paper, we studied the quantum effects generated by spherical partially transparent 
mirrors expanding with a constant acceleration in a D-dimensional flat spacetime. We 
consider a scalar massless field with an arbitrary parameter of non- minimal coupling £. 
The choice of this parameter does not affect the field equation and the expectation value of 
((p 2 ) Ten but results in different expressions for the stress-energy tensor. A partially trans- 
parent mirror is modeled by a 5-like potential in the field equation. Two special models 
were studied: Model A with a single mirror which expands with a constant acceleration 
so that its worldsheet is described by the equation R 2 — T 2 = b 2 , where b~ l is the value 
of the acceleration, and Model B with a couple of concentric expanding mirrors. 

We demonstrated that the leading terms of (ip 2 ) ien and (T^) rcn at J + has quite a simple 
form, ( |6.5| ) and ( |6.13| )-( |6.14| ), respectively. Both expressions infinitely grow at u — 0, the 
moment of the retarded time when the mirrors reach J + . The same divergence takes 
place at the moment v — of the advanced time when the mirrors start their motion 
from J~ . Both of these divergences are evidently connected with the adopted idealization 
of the problem: an infinite time of the accelerated motion and infinite sharpness of the 
mirrors. It should be also emphasized that since we used the Euclidean approach and 
obtained (0 2 ) vcn and (T^) rcn by using the Wick's rotation from their Euclidean values, the 
corresponding quantities in the physical spacetime are given for a special choice of state 
of the quantum field, namely the state which is invariant under time reflection T — > —T 
(for more details of the properties of such states, see p2| and [p3] ). 

It is interesting to note that in the model B for the conformally invariant case (that is, 
for a special choice of the coupling £ = {d— 1) / (4c?)), the stress-energy tensor identically 
vanishes both inside the inner mirror and outside the external mirror. For the conformal 



field in Model A this occurs everywhere. As indicated earlier (see [22] and [23]]), this is a 
direct consequence of the conformal invariance of the models. A more general discussion 
of the properties of the vacuum in Minkowski spacetime under conformal transformations 



can be found in [35] and 136] . 
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A Function F^ d \z 1 f3) and its properties 

In this Appendix we collect some useful formulae connected with properties of function 
F !yd \z,f3). This function is defined by the series 



It also allows the following representation (|32|, v.l, equation 5.2.11.15): 



o (1 — ze- 
Using this representation it is easy to see that 

i^(0,/3) = ^. (A.3) 

The function F^ d \z, (3) is defined for z < 1, where the series ( A.lj) is convergent. Near 
z = 1, it has a divergence. Since this divergence is connected with the behavior of the 
series at large L, one can estimate the leading divergence by putting (3 — 1 in (|A.1|) . 
Using the relation (||32||, v.l relation 5.2.11.1) 



(L + d-2)\ L 1 (d-3) 



—z 



1 - (1 - ^) d - 2 
~ (d - 2)z(l - z)^ 2 

we can conclude that the leading divergence of near z = 1 is of the form 



(A.4) 



In the main text, we also use the relation 



L!(d-2)! 

This relation allows one to show that 



1 l-f^^- (A,) 



lim 

,3->oo 



Using integral representation ( |A.2j ) one can verify that the following two relations are 
valid for functions F^(z,/3) 

dzFM(z,0) =-^-F^ d - 1 \z,(3-l), (A.8) 

(2 — 2 

dz F^ (z, (3) = — !— z d - x (z, 13) . (A.9) 

d 1 
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